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Abstract. Due to their high intrinsic brightness, caused by the disruption of the progenitor envelope by the 
shock-wave initiated at the bounce of the collapsing core, hydrogen-rich (type II) supernovae (SN) can be used 
as lighthouses to constrain distances in the Universe using variants of the Baade-Wesselink method. Based on a 
large set of CMFGEN models (Hillier & Miller 1998) covering the photospheric phase of type II SN, we study the 
various concepts entering one such technique, the Expanding Photosphere Method (EPM). 

We compute correction factors £ needed to approximate the synthetic Spectral Energy Distribution (SED) with 
that of a blackbody at temperature T. Our £, although similar, are systematically greater, by ~ 0.1, than the values 
obtained by Eastman et al. (1996) and translate into a systematic enhancement of 10-20% in EPM-distances. We 
find that line emission and absorption, not directly linked to color temperature variations, can considerably alter 
the synthetic magnitude: in particular, line-blanketing attributable to Fen and Tin is the principal cause for 
above-unity correction factors in the B and V bands in hydrogen-recombining models. 

Following the dominance of electron-scattering opacity in type II SN outflows, the blackbody SED arising at the 
thermalization depth is diluted, by a factor of approximately 0.2 to 0.4 for fully- or partially-ionized models, but 
rising to unity as hydrogen recombines for effective temperatures below 9,000 K. For a given effective temperature, 
models with a larger spatial scale, or lower density exponent, have a larger electron-scattering optical depth at the 
photosphere and consequently suffer enhanced dilution. We also find that when lines are present in the emergent 
spectrum, the photospheric radius in the corresponding wavelength range can be enhanced by a factor of 2-3 
compared to the case when only continuum opacity is considered. Lines can thus nullify the uniqueness of the 
photosphere radius and invalidate the Baade method at the heart of the EPM. Both the impact of line-blanketing 
on the SED and the photospheric radius at low T suggest that the EPM is best used at early times, when the 
outflow is fully ionized and line-opacity mostly confined to the UV range. 

We also investigate how reliably one can measure the photospheric velocity from P-Cygni line profiles. Contrary 
to the usually held belief, the velocity at maximum absorption in the P-Cygni trough of optically-thick lines can 
both overestimate or underestimate the photospheric velocity, with a magnitude that depends on the SN outflow 
density gradient and the optical thickness of the line. This stems from the wavelength-shift, toward line-center, 
of the location of maximum line-absorption for rays with larger impact parameters. This has implications for the 
measurement of expansion rates in SN outflows, especially at earlier times when only fewer, broader, and blue- 
shifted lines are present. This investigation should facilitate more reliable use of the EPM and the determination 
of distances in the Universe using type II SN. 
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1. Introduction 

In Dessart & Hillier (2005a, hereafter Paper I), we have 
presented first and general results obtained with the 
non-LTE model atmosphere code CMFGEN (Hillier & 
Miller 1998) for the quantitative spectroscopic analysis of 
photospheric-phase type II supernovae (SN). We demon- 
Send offprint requests to: Luc Dessart, 
e-mail: lucOas . ar izona . edu 



strated the ability of CMFGEN to reproduce, with great 
accuracy, the evolution of type II SN spectra. The only no- 
ticeable difficulty experienced by the code is in fitting Ha 
during the hydrogen-recombination phase, which is under- 
estimated both in absorption and emission by up to a fac- 
tor of 2-3 compared to observations. However, CMFGEN 
does a superb job for earlier times, both for the continuum 
and the lines, and from UV to near-IR wavelengths. This 
gives us confidence that CMFGEN can provide informa- 
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tion, at a very competitive level, on the properties of the 
SN ejecta and those of the progenitor. 

In this paper, we leave this aspect aside and focus on 
using CMFGEN to discuss current methods that use type 
II SN for distance determinations in the Universe. Being 
typically a factor of 10,000 times as bright as Cepheids, 
type II SN constitute an excellent tool to set distances 
over a larger volume of the Universe, while providing an 
alternative to type la SN and their assumed uniform prop- 
erties. 

Two methods are used to determine distances in 
the Universe based on type II SN: the Expanding 
Photosphere Method (EPM; Kirshner & Kwan 1974, 
Eastman & Kirshner 1989) and the Spectral-fitting 
Expanding Atmosphere Method (SEAM; Baron et al. 
1995, 2004). The EPM is based on the Baade method 
(Baade 1926), which requires a unique and well-defined 
radius for the expanding photosphere. As we shall see be- 
low, this assumption is only valid at early times when the 
outflow is fully ionized, i.e. when effects due to metal line- 
blanketing in the optical are small. While the EPM ap- 
proximates the Spectral Energy Distribution (SED) with 
a blackbody, adjusting for inconsistencies by means of cor- 
rection factors (Eastman ct al. 1996, hereafter E96; see 
Sect. EE} ) the SEAM uses synthetic spectra to directly fit 
the observed SED, and thus bypasses the use of correction 
factors. 

Distance determinations, based on the EPM, have 
shown different degrees of consistency with the Cepheid- 
based distance. The EPM works fine for the LMC (us- 
ing SN1987A, Eastman & Kirshner 1989) but predicts 
a distance to NGC 1637 50% smaller (using SN1999em, 
Leonard et al. 2002) than found using Cepheids (Leonard 
et al. 2003). On the contrary, every time the SEAM has 
been used, it has delivered the previously determined 
Cepheid distance. Why does the EPM work fine when 
applied to SN1987A and not for SN1999em? Indeed, the 
EPM should give an accurate distance, provided one uses 
the appropriate correction factors and photospheric veloc- 
ity. Thus of crucial importance to the EPM method is an 
understanding of the errors and biases associated with the 
derived correction factors and photospheric velocities. 

Due to its ease of use, and since it does not require flux 
calibrated spectra, it is desirable to investigate whether 
the EPM technique could be improved, and thus used re- 
liably for distance determinations. In this paper we de- 
vote our attention to the EPM and the various elements 
involved that may introduce errors into distance determi- 
nations: most importantly the correction factors, and the 
photospheric velocity measured from P-Cygni line profiles. 

This paper is structured as follows. In Sect. [21 we 
present the various aspects entering the EPM. In Sect. [31 
we study what correction factors we obtain in order to 
approximate the SED of type II SN with a blackbody. 
We also compare with equivalent values presented in E96. 
In Sect. ^ we investigate the physical origin and proper- 
ties of flux dilution by electron-scattering. We focus on 
one key quantity controlling the magnitude of such a di- 



lution, i.e. the ratio of thermalization to photospheric ra- 
dius. In Sect. we discuss how the model photospheric 
velocity compares with various measurements performed 
on P-Cygni line profiles in the optical range. Combining 
these various insights, we conclude in Sect. by present- 
ing the prerequisite for minimizing errors in the EPM and 
the Baade- Wesselink methods. 

2. Presentation of the EPM method 

Developed for spatially-resolved outbursts, the Baade 
method (1926) constrains the distance D to a novae based 
on measures of the rate of expansion i> p hot of the radiating 
photospheric surface of radius i? p hot, and of its apparent 
angular size 9 = R p hot/D, following the formula, 

D = v phot At/9, (1) 

where At represents the difference between the time of 
observation t and that of outburst to . Here we have made 
the key assumption of homologous expansion, which stip- 
ulates that mass shells, after a prompt acceleration at to, 
travel at constant velocity thereafter. We further neglect 
the original radius of the object, typically much smaller 
than corresponding values at times t of the order of days. 
When the outburst date to is not known, one uses a sam- 
ple of observations to fit the distribution of points along 
the line defined by t — (9/v p h ot )D + 1 , with slope D and 
ordinate crossing t . 

In principle, this method is directly applicable to a SN 
explosion suitably positioned within the Milky Way. These 
events are, however, too rare, and daily SN discoveries are 
instead permitted by searching over tens-of-thousands of 
galaxies. Such SN are too far to be resolved spatially, and 
thus, their angular size 9 must be obtained indirectly. 

For type II SN, the EPM estimates 9 by assuming that 
the outflow radiates as a blackbody at temperature T, 
yielding the relation, 

4^ hot ^(T) = 47ri? 2 /f rcd , (2) 

where all quantities (omitting D) are evaluated at fre- 
quency v. / dorcd is the de-reddened observed flux, i.e. 

Jdered = 1Q 0AA v fobs ^ w j th ^ thc cxtinction an( } jobs 

the observed flux. 

Equation (2), however, is ambiguous: the blackbody 
temperature T corresponds to the outflow layer where 
photons are thcrmalized (Rth), rather than the layer from 
where they escape (i? p hot), defined as the location where 
the total inward-integrated radial optical depth equals 
2/3. At early times of evolution, when the SN outflow 
is fully- or partially-ionized, i?th < fiphot because the 
electron-scattering opacity contributes significantly to the 
total opacity at the photosphere. Although this gray- 
opacity preserves the relative wavelength distribution of 
the Planck function ttB„(T), set at i?th, it introduces a 
global flux dilution. 

Schmutz et al. (1990) argued that, due to the finite 
extent of the photosphere, i?th covers a sizable range of 
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radii, increasing toward longer wavelengths, thereby asso- 
ciated with decreasing thermal temperatures and weaker 
dilution. When combined, these effects produce an emer- 
gent flux distribution characterized by a cooler blackbody 
temperature than that found at i? t h- Thus, to account for 
"sphericity" and dilution effects in Eq. @, we replace T 
by the color temperature T c , and i? p hot with £i? p hot, with 
explanations on how to obtain T c and £ given below. A 
more thorough discussion of the origin of the correction 
factor, £, is given in Sect.0] 

Taking these considerations into account , we can ex- 
press the observed flux as, 



/° bs - e 



9 2 lO" - 4 ^ 



nB v (T c ) . 



(3) 



The quantity £, often called dilution factor (e.g., E96; 
Hamuy et al. 2001, HOI), represents in fact a general cor- 
rection to be applied to the fitted blackbody distribution 
in order to reproduce the observed flux at v. 

By integrating each side of Eq. © over the transmis- 
sion function </>p it , of filter D, we now obtain an equation 
relating magnitudes, 



m D = -5 log(£) - 5 log(0) + A p + b T , 



where 



bn = -2.5 loe 



(4) 



(5) 



Here, we adopt the transmission functions 4>p tL , and zero- 
point constants C p from HOI. 

In practice, the color temperature T c is obtained by 
matching the slope of the corresponding blackbody to a 
set of magnitudes, usually a combination 5* taken from 
the optical and near-IR band passes UBVRIJHK . Not 
all band passes are equally suitable for use in the EPM: the 
fast-decreasing flux in the ?7-band makes the correspond- 
ing magnitude quickly too faint for good-quality observa- 
tion, so the J7-band is usually discarded; the i?-band also 
has problems due to its overlap with the Ha line, which 
shows significant variation in strength and width amongst 
type II SN. To allow a complete comparison between our 
results and those of E96, we sometimes quote in this work 
the results obtained for the near-IR band passes, although 
such magnitudes are rarely available. The various com- 
binations S usually adopted are thus, {B,V}, {B,V,I}, 
{VJ}, and {J,H,K}. 

The procedure is then, for a given set of photometric 
observations at time t, to apply a minimization of the 
quantity e defined as 



pes 



[m o + 51og[0£ s ]-A p -& o (T s )] S 



(6) 



For a given bandpass combination S and corresponding 
magnitudes, one obtains an ^-dependent color tempera- 
ture Tg: we thus write 9 as 6 s and £ as £g. 

For each date, one can thus determine the values of 8s 
and £s for all possible combinations of filters S as well as 



the corresponding value of u p hot from spectroscopic mea- 
surements. Grouping observations taken over a number of 
nights, one can then plot the quantity, 



t = D 



0s_ 

Wphot 



+ *0 



(7) 



Fitting a line to the resulting distribution of points deliv- 
ers the distance D (the slope) and the time of explosion 
(the ordinate crossing). One can compare the results for 
different choices of band passes S. 

A fundamental part in the EPM is to quantify, in 
Eq. ©, the correction factor £s associated with the black- 
body color temperature Tg, for each bandpass set S. The 
procedure is to recast Eq. in terms of absolute magni- 
tude and distance modulus, 



M 



-51og£s - 5 log 



-Rphot 

10pc 



(8) 



Then, based on a large set of model atmosphere cal- 
culations providing measures of -Rphot and Mp, one can 
minimize, for the different bandpass choices 5, the quan- 
tity e given by, 



E 

pes 



Af p + 51og[£s] +5 log 



-Rphot \ 

10pc ) 



b e (T s 



(9) 



E96 provides analytical formulations to describe how 
£s changes with Tg, subsequently revised by HOI to in- 
clude a more suitable set of filter transmission functions 
and zero-points. Such formulations have been used by ob- 
servers to obtain distance estimates to type II SN (see, 
e.g., Schmidt et al. 1994ab; E96; HOI; L02). In Sect. 
based on CMFGEN models, we compute similar relations 
to compare and gauge the uncertainties, since these trans- 
late directly into the EPM-distance. In Sect. 01 we make a 
digression on the physical origin of flux dilution, its prop- 
erties for a wide variety of type II SN and its correspon- 
dence with correction factors. In Sect.[Sl we then discuss 
the procedure to infer the photospheric velocity, a quan- 
tity whose uncertainty weighs on the EPM just as much 
as correction factors do. 

3. Correction factors 

The EPM-distance, obtained by minimizing e (Eq. EJ), re- 
quires knowledge of the dependence of correction factors 
£s on blackbody color temperature Tg, whose analytical 
description, at present, is only given in the work of E96. 
Here, following E96 and the casting of e given in Eq. (JSJ, 
we compute the variation of £s with Tg based on a large 
number of CMFGEN models that cover the type II SN 
evolution over the first few weeks after core-collapse; then, 
we compare with the results of E96. 

In Figs. □ FJ and we plot various combinations of 
the resulting quantities that minimize this last form of e. 
Figure H shows the variation of the correction factors £s 
versus the blackbody color temperature, Ts, obtained for 
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each of the four bandpass combinations, {B, V}, {B, V, I}, 
{V,/}, and {J,H,K}. Figure |21 shows the variation of 
Ts/T e R versus the effective temperature T c g, defined here 
as the gas (electron) temperature at continuum optical 
depth two-third (integrated inwards from the outer grid 
radius). This location corresponds to the photospheric ra- 
dius -Rphot- Finally, Fig. [3] shows the variation of £s with 
photospheric density p p h- We refer the reader to E96 for a 
discussion of the properties shown in Fig. [3 and 01 (which 
are similar to theirs), and here focus our attention to the 
relation £s-Ts illustrated in Fig. Q 
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Fig. 1. Correction factors £s versus the corresponding 
blackbody color temperature Ts, obtained simultaneously 
by minimizing e (Eq. . We show solutions for four differ- 
ent bandpass combinations, {B, V}, {B, V, I}, {V, I}, and 
{J,H,K}. In each panel, we overplot a first- or second- 
order polynomial of the form £g t = J2i a i(10 4 K/Ts)\ with 
coefficients obtained from a best-fit to these data points 
(solid line; see Tabled or adopted from HOI (dashed line). 
For ease of comparison, a unique ordinate range is used 
for all plots. 

Let us focus on the top-left panel of Fig. ^ showing 
the variation of £gy as a function of blackbody color tem- 
perature Tbv- The distribution of data points follows two 
distinct patterns: above ca. 9,000 K, £gy is roughly inde- 
pendent of Tbv and of order 0.45, but below ca. 9,000 K, 
it increases for lower values of Tbv, culminating at 1.5 for 
Tbv ~ 4000 K. For a fixed temperature Tbv, the scatter 
is uniform and of order 0.1, corresponding to decreasing 



Fig. 2. Solutions obtained from Eq. (0, this time show- 
ing the color temperature in each bandpass combination 
S normalized by, and plotted as a function of, T e gf. The 
band passes are arranged as in Fig. 1, and for ease of com- 
parison, a unique ordinate range is used for all plots. 

uncertainties towards lower temperature (from 20 down 
to 10%). But over a range of temperature, this scatter in- 
creases, at low Tbv, by a factor 2-3, following the sharp 
rise of £,bv- Thus, adopting a unique correction factor at 
a given temperature can introduce non-negligible errors, 
either through the neglect of the intrinsic scatter of £bv, 
or through an inaccuracy in the temperature Tbv, made 
worse for smaller values of Tbv- In addition, we overplot 
a polynomial with coefficients obtained from a best-fit to 
our distribution of data points (solid line; see TableQ]), or 
extracted from Table 14 in H01 (dashed line). Despite a 
qualitative agreement, our best-fit polynomial lies system- 
atically above HOl's or E96's, with a maximum difference 
of ca. 0.2 (or 50%) at T BV ~ 8,000 K. 

Turning to the other panels in Fig. ^ we see that £sy / 
and £yj follow a similar trend, although £yj remains be- 
low unity at low temperatures Tyj. For £,jhk, the trend is 
altogether different, correction factors decreasing at lower 
temperatures, with a scatter reduced compared to the pre- 
vious three cases. Note that the {J,H,K} bandpass is 
poorly sensitive to color temperature in the 10,000 K tem- 
perature range. Here again, the overplotted polynomials 
show that our data points lie, on average, above those 
obtained by H01/E96. 

The reason for the discrepancy is unclear. It might be 
related to the different approach E96 adopts to handle 
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Fig. 3. Same as Fig.^showing the correction factors as a 
function of photospheric density, again for each bandpass 
combination S. For ease of comparison, a unique ordinate 
range is used for all plots. 

the relativistic terms. Also, E96 only solved the non-LTE 
problem for a few species, obtaining a consistent distri- 
bution between absorptive and scattering opacity contri- 
butions, while for the rest, usually metals, line opacity 
was taken as pure scattering. It is not clear that this 
can cause a sizable relative shift of the thermalization 
and photospheric radii, controlled by hydrogen bound- 
free and electron-scattering opacities (both handled in 
a similar way in E96 and here). Finally, E96 handles 
metal line opacity with the "expansion opacity" formal- 
ism of Eastman & Pinto (1993), a convenient but approxi- 
mate approach compared to the more consistent treatment 
adopted in CMFGEN, where all line opacity is included 
systematically in a "brutal" but accurate way. The only 
other obvious difference in our approach is the use of a very 
large model atom (manageable through the use of Super- 
Levels) , to allow a detailed and complete census of opacity 
sources, as well as the non-LTE treatment of all species. 
These resulting differences are not academic because our 
systematically larger correction factors translate directly 
into EPM distances that are, on average, 10-20% larger 
than predicted by H01/E96 (Dessart & Hillier 2005c). 

Let us now leave the discussion on correction factors 
and investigate the adequacy of blackbody color tempera- 
tures obtained from the minimization of e given in Eq. JSJl . 
In Fig. we plot the blackbody SED (solid, featureless 
line) that best fit the SED of a CMFGEN model (T off = 




X (A) 




X (A) 




X (A) 




Fig. 4. Comparison between the model with 
T cff =9,200 K that fits well the observations of SN1999em 
on the 5th of November 1999 (Sect. 4) with black- 
body spectra whose temperature is obtained by fitting 
the model magnitude in different band passes {B,V}, 
{B, V, I}, {V, I}, and {J, H, K}, shown in this order from 
the top to the bottom panels. The corresponding color 
temperatures and correction factors are: 12,945 K and 
0.443 (BV), 11,951K and 0.494 (BVI), 11,278K and 
0.526 (VI), 7,174K and 0.843 (JHK). The dotted line 
shows the spectral energy distribution of the model when 
only continuum processes are included. 
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Table 1. Coefficients of the polynomial fits (shown as 
solid lines in Fig. ^| to the distribution of correction fac- 
tors £s versus blackbody color temperature for the four 
bandpass combinations. As in H01/E96, we use a second- 
order polynomial of the form £fi t = a i (10 4 K/Ts) 1 . 





{B,V} 


{B,V,I} 


{v,i} 


{J,H,K} 


a 


0.47188 


0.63241 


0.81662 


0.10786 


ai 


-0.25399 


-0.38375 


-0.62896 


1.12374 


a 2 


0.32630 


0.28425 


0.33852 


0.00000 




X (A) 



9,200 K; solid line) over bandpass combinations {B,V}, 
{B, V, I}, {V, I}, and {J, H, K}. Blackbody color temper- 
atures (correction factors) derived for each set are, in the 
same order, 12,945K (0.443), 11,951 K (0.494), 11,278K 
(0.526), and 7,174 K (0.843), thus above and below T cff 
(Fig. [2Jl - The fit quality is good in all cases, although the 
blackbody SED shows a sizable vertical shift, attributable 
to the presence of lines, above the synthetic continuum 
SED (dotted line). Indeed, Mp, in Eq. accounts indis- 
criminately for the total flux in the corresponding band- 
pass, be it continuum, line, or a combination of both. If, 
however, one uses, in Eq. (0, magnitudes computed on 
the synthetic continuum SED (Fig. [SJ, the correspond- 
ing blackbody SED then fits accurately the synthetic 
continuum SED. The color temperatures are then lower 
and equal, in the same sequence, to 11,599K, 10,957K, 
10,316 K and 6,565 K, with associated correction factors 
5-10% bigger. 

The influence of lines on the optical SED strengthens 
as the outflow temperature further decreases. A model 
with T c ff =5,950 K illustrates more dramatically this sit- 
uation in Fig. [fj] while the blackbody SED fits well the 
observations in the {J, H, K} set, the agreement becomes 
poor for the sets {V, 1} and {B, V, I}, and unsatisfactory 
for the set {B, V}. Here, the severe line-blanketing due to 
Fen and Tin (Paper I) causes a major flux deficit around 
4200A; the resulting blackbody color temperature, Tbv = 
4,545 K, is too low to represent well the optical SED (note 
the absence of a clearly defined continuum) and correc- 
tion factors, which adjust according to constraints from 
Eq. (O, rise dramatically, as do the associated uncertain- 
ties. Despite its weak influence on the effective tempera- 
ture T ff of the outflow, an enhanced Ti abundance could 
thus reduce the B-band flux and the inferred Tbv , mimick- 
ing a distance dilution. Retaining observations sampling 
the early evolution of type II SN, i.e. prior to the ap- 
pearance of metal lines in the optical spectral range, will 
therefore reduce the uncertainty on distances inferred with 
the EPM. 

Within these limitations, the EPM should however de- 
liver a distance to the object with an accuracy set by the 
level of agreement between the synthetic and the observed 
fluxes (or magnitudes), since Eq. @ provides the neces- 
sary correction to apply to a blackbody SED to retrieve, 
within numerical accuracy, the synthetic magnitudes for 



{BVI} 



4000 6000 
A (A) 




4000 6000 
\ (A) 




Fig. 5. Same as previous figure but this time taking the 
continuum energy distribution. The corresponding color 
temperatures and dilution factors are now: 11,599 K and 
0.486 (BV), 10,957K and 0.529 (BVI), 10,316K and 
0.567 (VI) and 6,565 K and 0.895 (JHK). Note that both 
dilution factors and blackbody color temperatures are dif- 
ferent by about 10% from the previous case where the 
synthetic spectrum contained lines. 



any bandpass combination. Overall, the above argument 
suggests that the EPM will be more accurate when com- 
bined with a detailed multi-wavelength modeling of, rather 
than a degenerate 2- or 3-point magnitude fit to the ob- 
served SED, thereby lifting possible degeneracies. A fur- 
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Fig. 6. Blackbody fits in the four band passes to the 
model that reproduces well the observations of SN1999em 
taken on the 22nd of November 1999 (see Dessart & Hillier 
2005c), i.e. at ca. one month past core-collapse. The cor- 
responding color temperatures and dilution factors are 
now: 4,577 K and 1.471 (BV), 5,346 K and 0.914 (BVI), 
6,180K and 0.7 (VI), 5,635K and 0.790 (JHK). Note 
how inadequate is the fit to the BV set, resulting in a 
correction factor greater than unity. Conceptually, correc- 
tion factors are not dilution factors since they correct for 
a number of inadequacies of the EPM, presently the cor- 
rupting effect of lines. 



ther advantage of multi-wavelength modeling is that the 
influence of lines can be more accurately gauged and hence 
allowed for. 

4. Dilution factors 

We now investigate the elements that contribute to 
the flux dilution caused by the dominance of electron- 
scattering opacity. In Sect. 14.11 we discuss its origin in 
more details. A key quantity that controls the level of 
dilution is the spatial shift between the thermalization ra- 
dius (Rth) and the photospheric radius (i? p hot)- We show 
in Sect. l4.2l how these different radii vary with wavelength, 
and provide an additional discussion on the effect of line- 
blanketing. In Sect. 14.31 we extend the discussion to a large 
set of CMFGEN models, drawing statistical properties of 
Rth/Rphat over a wide parameter space, i.e. corresponding 
to different density exponents, outflow ionization, outflow 
velocity etc... 

4.1. Dilution factors and their origin 

The origin of dilution factors has been discussed several 
times in the literature (E96; Hershkowitz & Wagoner 1987; 
Chilikuri & Wagoner 1988; Schmutz et al. 1990). In order 
to clarify some concepts we repeat some of this discussion. 
First, consider a plane parallel scattering atmosphere, and 
define A as the ratio of absorptive opacity K a b s to the total 
opacity K to t (i.e., A = Kabs/«tat)- If we assume that A is 
constant with depth, and that the Planck function varies 
linearly with r (i.e, B(t) = a + br) it can be shown that 
the mean intensity J is given (approximately) by 



J(t) =a + br- 



b — \/3a 

7T- 



exp[ 



3A- 



(10) 



(e.g., Mihalas 1978). We define the thermalization optical 
depth, A, by 1 



A = 0.67/VA 



(11) 



Using this definition, and the previous expression it is ap- 
parent, for arbitrary a and b, that J will depart signif- 
icantly from B when the optical depth is less than the 
thermalization depth, i.e. r < A. 

It can be further shown that the astrophysical flux F 
is given, approximately, by 



F 



V31 + VA 



B(t = 1/V3A) 



(12) 



Thus the observed spectrum is that of a blackbody at 
r = 0.87A but diluted by a factor 



(13) 



1 Mihalas defines the thermalization depth with unity in- 
stead of 0.67, while Schmutz et al. (1990) use 0.58 (= 1/V3). 
We adopted 2/3 because we use that optical depth to define 
the photospheric radius. These slightly different definitions are 
of no practical consequence. 
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For a strongly scattering atmosphere (i.e, A << 1), the 
dilution (= £ 2 ) is approximately 2\/A. 

Unfortunately, the situation with a SN is much more 
complicated. First, A is depth dependent, since the scat- 
tering opacity (due to free electrons) scales linearly with 
density while the absorptive opacity, due to bound-free 
and free-free processes, scales typically with the square of 
the density. However, we can define an approximate ther- 
malization optical depth scale by 

/oo 
VXKtotdr. (14) 

The thermalization layer is located on this optical depth 
scale at a value of 2/3. Since A varies with wavelength, 
r t h also does. Consequently the observed spectrum is not 
that of a pure blackbody. 

Second, SN atmospheres have a finite extent. Thus, 
the photosphcric radius i? p hot (defined, using the plane- 
parallel results, as occurring at r = 2/3) is larger than 
the radius of the SN at the thermalization depth (i?th)- 
This introduces an additional dilution of the radiation — 
roughly as (i?th/-R P hot) 2 - 

Third, the finite extent of the envelope also introduces 
an additional correction due to emission from the overly- 
ing layers. 

Given these complications it makes sense to write, as 
done in Eq. (3), 

/° bs = e 2 6 2 10- 4A - ttB v (T c ) (15) 

where £ represents a factor that corrects for the measured 
blackbody flux to give the observed flux. Note that with 
this approach, the precise definition of photosphcric ra- 
dius is unimportant — what is important is that we com- 
pute the £ for whatever definition we use. As discussed in 
Sect. |3 the correction factors entering the EPM are not 
simply accounting for the effects just discussed, but also 
for the significant impact that line-blanketing has on the 
SED. 

In the next section we illustrate how the thermalization 
and photosphere radii vary with wavelength. In terms of 
the thermalization radius we find that the total dilution 
(£ 2 ) scales as 

(i?th/-Rphot)" +1 

where n is the density exponent in the SN, and we have 
assumed that the ratio of absorption to scattering opacity 
is constant. Thus, a small change in thermalization radius 
can produce a substantial change in the dilution factor. 

4.2. Wavelength variation of and i? p hot in realistic 
SN models 

We have selected three models covering the early-, 
intermediate- and late-stage evolution of a type II SN dur- 
ing its photospheric phase. In the top panel of Fig. [7| we 
plot the wavelength variation of the radius of optical depth 
two-third for a model with T c g — 13,6 10 K, normalized to 
the base radius Rq (also noted R*) of the model grid. 
Three curves are illustrated: 



1. R(t = 2/3) when all processes, such as bound-bound, 
bound-free, free-free, and electron-scattering are in- 
cluded (black curve). 

2. R(t = 2/3) when all processes, except bound-bound 
transitions, are included (blue) . This defines the classic 
photosphere. In most regions, it is set by the electron 
scattering opacity. 

3. i?(r t h = 2/3) where r t h is computed using Eqn. ifbljl 
and we assume that the scattering opacity is due to 
electrons, while the absorptive opacity is due to bound- 
free and free- free processes only (red). 

Several elements are apparent in this plot. First, the 
curves depart from the horizontal by various degrees, i.e. 
the variation with wavelength depends on the processes 
included. The thermalization radius, although generally 
increasing with wavelength, shows a sawtooth behavior. 
The main opacity process that links photons to the ther- 
mal pool is the photo-ionization of hydrogen, the sud- 
den jumps in thermalization radius at 3649A (Balmer 
jump) and 821lA (Paschen jump) resulting from photo- 
ionization from the n = 2 and n = 3 levels of hydrogen. 
The thermalization radius varies between 1.1 and 1.25 Rq 
where Rq is the base radius of the model. While this 
represents a modest radius variation, the corresponding 
electron-temperature decreases from 25,000 K to 17,000 K 
between these two heights so that the associated Planckian 
photon distribution at 1000 A and 25000 A at the thermal- 
ization depth are in fact substantially different. The sec- 
ond curve up is remarkably wavelength independent, fixed 
at 1.4 Rq. It corresponds to the photosphere, for which one 
usually assumes only continuum processes. 

We show in the bottom plot the rectified flux over the 
same wavelength range as that of the top panel. Note how 
the flux minima and maxima are "echoed" in the top panel 
with large increases in radius. This is particularly true 
in the UV where metal line-blanketing around 2000A is 
strong. Thus, strictly speaking, the spatial size of the pho- 
tosphere is very dependent on wavelength, but this vari- 
ation has different magnitudes depending on the selected 
opacity sources and the wavelength region considered. 

In Figs. 0D and we show a similar set of plots but 
now for models of decreasing effective temperature (out- 
flow ionization), corresponding to a later stage of SN evo- 
lution where hydrogen recombines at or above the photo- 
sphere. One key feature shown in these plots is that the 
thermalization radius is now much closer to the photo- 
spheric radius: the reduced density of free-electrons favors 
the contribution of bound-free and free-free processes to 
the total optical depth at the photosphere. 

A second effect clearly appears in these latter two fig- 
ures, namely the presence of lines, first predominantly in 
the UV (Fig.[8J), and then also in the optical range (Fig.^J). 
The effects of metal lines in the UV increases the photo- 
sphere radius from ca. 1.3 Rq up to ca. 2-4 Rq. In the opti- 
cal range, the B and to a lesser extent the V bands become 
more and more contaminated by lines so that the effective 
photospheric radius is 10-15% bigger than predicted when 
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Fig. 7. Top panel: wavelength variation of the radius at 
optical-depth two-third, including from the bottom to the 
top curves, only bound-free and free- free processes (red), 
bound-free plus free-free plus electron-scattering contri- 
butions (blue), and all continuum and line processes to- 
gether (black). The effective temperature of the model is 
T e ff =13,610K. Bottom panel: corresponding wavelength 
variation of the rectified flux. Also overplotted are the 
transmission curves of the B, V, I, J H, and K bands 
from Bessell (1990, 1988) and HOI. The rectified flux was 
obtained by dividing the model spectrum by the model 
continuum spectrum. 



accounting only for continuum opacity. So, the first ef- 
fect of lines at later times is to spread the photosphere 
radius over a range of values, making it non- unique and 
therefore less adequate for the application of the Baade- 
Wesselink method. But as the outflow cools down, lines 
become ever stronger and it becomes more and more du- 
bious to approximate the SED with a blackbody distribu- 
tion of whatever temperature (EPM). In the I band, which 
at early times was completely devoid of lines, the Can line 
around 8800A, causes a huge variation in the (true) pho- 
tosphere radius of cooler models. It is again tempting to 
advise not to use the I band at later times for the EPM. 

Let us illustrate further the strong impact lines can 
have on modifying the photospheric radius. We show in 
Fig. contours of total optical depth two-third in the 
(z,p) plane (in units of Rq) at selected positions within 
the Ha line profile, placing the observer at infinity to the 
left [p = 0, z = —00). The black half-circle corresponds to 
the base of the grid at Rq. Here, the optical depth includes 
all processes, but over the small wavelength region con- 
sidered, only the contribution from the line optical depth 
varies. Each contour corresponds to a given wavelength 



Fig. 8. Same as Fig. [7| for a cooler model (intermediate- 
stage in the photospheric-phase evolution), with 
T e ff =8,035 K. The maximum ordinate in the top 
panel is 3, rather than 2.3 as in Fig.0 




10000 15000 
A (A) 



Fig. 9. Same as Fig. for a much cooler model (late- 
stage in the photospheric evolution), with T c ft =6,460 K. 
The maximum ordinate in the top panel is 3, rather than 
2.3 as in Fig. EJ 



displacement from line center, shown as a vertical line of 
the same color in the right panel. Throughout the profile, 
the minimum optical depth is fixed by continuum pro- 
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cesses, which, as said before, are essentially constant over 
the line profile. This corresponds to the curve that is clos- 
est to the inner half-circle, shown in red. Any departure 
from that curve is due to line-optical depth effects. 

In outflows characterized by a Hubble velocity, the iso- 
velocity contours (i.e, contours on which the projected line 
of sight velocity is constant) coincide with z=constant 
(i.e., planes perpendicular to our line of sight). Thus, if 
the line optical depth was only dependent on velocity, the 
iso-velocity and iso-r contours would overlap. However, 
the line optical depth decreases with the density, which 
drops as a high power of the radius in SN outflows. Thus, 
an iso-velocity (or constant-z) surface cuts iso-density con- 
tours for increasing values of p. Beyond a given maximum 
impact parameter, the line optical depth becomes negligi- 
ble and the total optical depth takes its minimum value. 
For sufficiently high values of p, there is no surface of total 
optical depth unity. The last such surface defines the limb 
of the star. 

At blue- wavelengths from line center, the iso-r contour 
is strongly peaked along p = 0. At wavelengths closer to 
line center, the iso-r surfaces cover a bigger range of im- 
pact parameters, going up to a value of 2.3 i?o which cor- 
responds to the wavelength for maximum flux (line center 
region). This means that the radiating surface at 6562. 79A 
is 70% bigger than that just redward of Ha where no 
lines are present. We thus see that the radiating surface 
has a wavelength-variable radius once lines are accounted 
for, invalidating further the assumption of constant pho- 
tospheric radius and blackbody energy distribution. These 
facts support the idea that the Baade-Wesselink method 
will work best during the early-stage of the photospheric 
phase evolution of type II SN, when they show a rela- 
tively featureless spectrum and a more uniform dilution 
from electron scattering. 

4.3. Statistical properties of Rth/Rphot in realistic SN 
models 

In the previous section, using a set of three CMFGEN 
models of distinct outflow ionization, we have discussed 
the wavelength variation of the location of selected radii, 
corresponding to thermalization, continuum and total op- 
tical depth two-third. If lines are ignored, the photospheric 
radius is essentially unchanged throughout the spectrum. 
If lines are accounted for, as they should be, the result- 
ing "true" photospheric radius can vary widely over small 
wavelength ranges. For the following discussion we make 
the important simplification that the photosphere is only 
set by continuum opacity processes. This will be inaccu- 
rate at intermediate- and late-times but reasonable during 
the early-phase of photospheric evolution of type II SN 
(i.e. until metal- lines appear in the optical range). Below 
we investigate how the quantity i?th/-Rphot varies over a 
wide range of outflow properties, controlled, e.g., by out- 
flow ionization, the density distribution, and the spatial 
scale. Using the modest variation of i? t h with wavelength, 



we simplify the discussion by building optical and near-IR 
averages (i.e. over the sets {B,V,I} and {J, H, K}) for 

-Rth/-Rphot- 

4.3.1. Dependence on T c g and p p hot 

In Fig. 1111 we first illustrate the variation of -Rth/-Rphot 
with T c ff and p p hot, using averages over the two sets 
{B,V,I} and {J, H, K}. Let us focus on the {B, V, 1} set 
and discuss the results for the n = 10 models shown as 
green dots (we leave the discussion for the dependence on 
n until the next section). 

For this sub-set, there is a very modest scatter of 
Rth/Rphot in the temperature range 8,000 to 17,000 K, 
by ±0.03 around ca. 0.68. At the low temperature end, 
this ratio rises towards unity. This is easily understood 
by recalling the previous section. Above approximately 
8,000 K, the outflow temperature and ionization is such 
that most electrons are unbound and thus contribute a 
significant and more or less fixed fraction to the total 
continuum optical-depth at the photosphere. At approxi- 
mately 8,000 K, hydrogen suddenly recombines in the pho- 
tosphere region and the density of free-electrons drops. For 
lower photospheric temperatures, the electron-scattering 
contribution to the total continuum optical-depth at the 
photosphere decreases in favor of bound-free and free-free 
opacity sources, so that the corresponding dilution fac- 
tors increase towards unity. This limiting case is reached 
when electron-scattering no longer contributes to the op- 
tical depth at the photosphere, i.e. when no free-electrons 
are present. In this case, the SED at the photosphere suf- 
fers no dilution; indeed, it is a fundamental property that 
-Rth/-Rphot be less than unity. 

When Rth/Rphot is plotted as a function of photo- 
spheric density, the models with n = 10 (green-dots) also 
show a tight correlation, falling on a positive-gradient 
line with increasing density. As the density increases, 
-Rth/-Rphot increases towards unity. This occurs since a 
higher and higher photospheric density is required to al- 
low bound-free and free-free processes to compensate for 
the reduced contribution of electron-scattering at the pho- 
tosphere (due to recombination of hydrogen). Moving to 
the lower panel and the corresponding illustrations for the 
set {J, H, K}, we find that i?th/^phot shows a very similar 
behavior, with values being merely shifted upwards by up 
to ca. 0.1. 

The tight correlation of -Rth/-Rphot with photospheric 
density and effective temperature is not clearly interpreted 
since these quantities may not be completely indepen- 
dent. To disentangle potential combined effects, we show 
in Fig. El a contour plot for i?th/-Rphot averaged over the 
{B,V,I} set, and limited still to models with n = 10 
(green dots in Fig. ^] note also that similar findings 
would be obtained for the near-IR set). As before, we find 
that in the fully ionized models there is a pronounced 
degeneracy of i?th/-Rphot with temperature, as revealed 
by near- horizontal contours. Points with i?th/-Rphot 1 
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Fig. 10. Left: Optical-depth-two-third contours for Ha, shown in the (z,p) plane, for an observer located at p = 
and z = -co. Each contour corresponds to a given Doppler-shift velocity from the rest wavelength of Ao = 6562. 79A. 
Right: Emergent Ha profile shown in velocity space, i.e. scaled wavelength displacements from the rest-wavelength 
(A/Ao — l)c. To assist with the visualization, each contour on the left corresponds to a vertical line of the same color 
in the right panel that gives off the selected velocity away from line-center. The red curve is for the continuum which, 
for practical purposes, is constant across the line [color]. 



occur in the parameter space where the T e g is low and 
Pphot is high, i.e. higher i?th/-Rphot are found as one moves 
from the lower right corner (low p p hot and high T e s) to 
the upper left corner (high p p hot and low T c «). As dis- 
cussed above, as one moves from fully ionized models to 
cooler ones where hydrogen recombines in the outflow, 
the electron-scattering optical depth at the photosphere 
is significantly lower due to the reduced density of free- 
electrons. For reduced photospheric temperatures, higher 
photospheric densities are required to compensate for the 
reduced contribution of electron-scattering to the contin- 
uum optical depth, explaining the rise of Rth/ R v hot from 
the bottom-right corner to the upper- left corner of Fig.lT2*l 



4.3.2. Dependence on secondary parameters 

Let us now broaden the discussion of the previous section 
by investigating the dependence of -ffth/-Rphot on model 
parameters at large, and not just the density and temper- 
ature at the photosphere. First, let us go back to Fig. 1111 
and consider all points. The scatter at a given temper- 
ature is then about 3-4 times as large as in the case 
where a single color is considered. When the color-coding 
is employed, models with a given density exponent are 
grouped together and the resulting scatter of points is 
quite small. The effect is easily understood if one considers 
the ideas discussed in Sect. 4.5 of Paper I. Thermalization 
processes are density-square dependent, while electron- 
scattering is linear in the density. Lowering (increasing) 
n makes the density in the outer layers larger (smaller) 
and thus leads to a lower (larger) density at the pho- 
tosphere: the opacity due to thermalization processes is 
lowered (increased) and i?th/-Rphot is reduced (enhanced). 
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Fig. 11. Plot of i?th/-Rphot averaged over the two sets 
{B, V, 1} and {J, H, K} versus effective temperature (left) 
and photospheric density (right) for 142 CMFGEN models 
covering a range of outflow ionization pertinent for model- 
ing the photospheric evolution of type II SN. To partially 
explain the scatter of points, we use a color-coding that 
groups models that have the same density exponent (see 
Sect. I4.3.1I4.3~2*|> . Note that a number of models which 
do differ in spectroscopic appearance have comparable 
-Rth/-Rphot, leading to numerous overlap of data points. 
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In general, a lower n favors the contribution of electron- 
scattering in fully ionized models and therefore results in 
low values of -Rth/-Rphot- ln c °ol models, the density expo- 
nent has a less pronounced importance since the outflow 
optical depth is essentially set by the optical depth of the 
hydrogen-recombination front. This typically extends over 
a spatially confined region, while the outer regions con- 
tribute a negligible electron-scattering optical depth due 
to the low density of free-electrons. 

In the right-hand side of Fig. ^2 a number of mod- 
els fall below the curved line followed by most models. 
These correspond to models that have a different scale, 
i.e. smaller luminosity and radius but similar outflow ion- 
ization (same T c ff). For example, setting the base radius 
at Rq and 10 Ro, the base density at po and po/10 and the 
luminosity as Lq and 100 L$ gives two models with iden- 
tical photospheric velocity and effective temperature and 
identical electron-scattering optical depth at the grid base. 
However, the model with a larger radius has i?th/-Rphot 
smaller by ca. 10% (enhanced flux dilution) than that of 
the model with the smaller radius. This is explained in the 
same terms as above for the models with different n. For 
more compact objects, the reduced scale of the outflow 
reduces the contribution from electron scattering (which 
varies linearly with the spatial scale) and leads to a higher 
density at the photosphere. Also, we find no noticeable 
metallicity effect on the resulting i?th/-Rphot- This is easily 
understood — in type II SN the thermalization processes 
stem from the photoionization of hydrogen atoms which is 
unrelated to the metal opacity. Moreover, modulating the 
model expansion velocity does not change the behavior of 

Rth/Rphot- 

Finally, in Fig. we duplicate Fig. fR] by including 
all the models shown in Fig. fT^ The scatter is signifi- 
cantly enhanced, although the general trend follows the 
same rule as observed with models having n — 10 only. 
So, while small variations from model to model can result 
from changes, mostly, on the density exponent and the 
spatial scale, the general behavior of i?th/-R P hot is always 
towards higher, i.e. closer-to-unity, values for cooler tem- 
peratures and concomitant higher photospheric densities. 

5. Inference of the photospheric velocity 

To infer the distance to a type II SN, the EPM requires 
the evaluation of the photospheric velocity at each ob- 
servation date. The standard procedure is to adopt the 
velocity at maximum absorption u a b s in P-Cygni profile 
troughs, favoring the more optically-thin lines (L02). This 
measurement has an inherent uncertainty, particularly at 
early times when profile troughs show a weak curvature; 
besides, it is unclear how suitable this location is since line 
profile formation and shape result from complex line and 
continuum optical-depth effects (Sect. 5, Paper I). In this 
section, we address the reliability of the above procedure 
by comparing synthetic line-profile velocity measurements 
with the corresponding model photospheric velocity. 
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Fig. 12. Contour plot of < i?th/-R P hot >{b,v.i} in the 
T e g — pphot plane for models with n — 10, shown as green 
dots in Fig. fT^ Iso-contours are this time shown for values 
of 0.80, 0.85, 0.90, and 0.95. 
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Fig. 13. Contour plot of < i?th/i? P hot >{b,vj} in the 
Tcft — Pphot plane. For visibility purposes, we only show 
iso-contours for ^-values of 0.7, 0.8, 0.9, and 0.95. 

In Fig. 1141 using a sample of ca. one hundred 
photospheric-phase type II SN models, we show the ve- 
locity measurement at maximum absorption in synthetic 
line profiles of Ha, H/3, H7, and Fell 5169A, normalized to 
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Fig. 14. Plot based on measurements performed on syn- 
thetic spectra of fabs/^phot as a function of i>phot for Ha, 
H/3, H7, and the Fen 5169A feature. i> a b s corresponds to 
the velocity at maximum absorption in the corresponding 
P-Cygni line profile. v p hot is the model velocity at which 
the total continuum optical depth is 2/3. Note that there 
are fewer data points for the Fe n feature since it is only 
observable in relatively cooler outflows, which statistically 
possess a lower expansion velocity [color] . 




Fig. 15. Plot of line (thick lines; dashed-dotted: Ha; solid: 
H/3; dashed: H7) and continuum (thin line; solid: electron- 
scattering; dashed: continuum at 1000A; dashed-dotted: 
continuum at 3000A; dotted: continuum at 5000A) op- 
tical depth as a function of velocity for the model dis- 
cussed in Sect. 03 The photosphere is located at a velocity 
of 8750 km s" 1 , at which the optical depth in hydrogen 
Balmer lines is orders of magnitude larger than that in 
the continuum. 



and as a function of the model photospheric velocity v p hot 
- statistically, models that have a higher v p hot also have a 
higher T c g and correspond to times closer to the explosion 
date. If the velocity location at maximum absorption in 
the synthetic line profile corresponded accurately with the 
model photospheric velocity, all data points would gather 
along the v a bs/« P hot = 1 line. Instead, we see a moderate 
scatter, between 0.8 and 1.4, enforced by the strong con- 
finement within the photospheric layer of line and contin- 
uum formation (Sect. 5, Paper I). Nonetheless, this scat- 
ter, a function of the synthetic line chosen, is large from 
the perspective of the EPM (Sect. with points both 
above and below this unity-line, implying a potential over- 
estimate or underestimate of the photospheric velocity. For 
Ha, most data points lie above the unity-line and thus 
adopting u a b s for that line generally overestimates v p hot, 
by up to 40%, but only marginally for high-v p hot /high-T e g 
models. For H/3 and H7, data points show a similar scat- 
ter as Ha, but shifted down by ca. 20%, matching v p h t 
by ±15%. But for the Fen 5169A line, the velocity mea- 
surement i> a bs matches w p hot to within 5-10%. 

This suggests that, when applying the EPM or the 
SEAM, setting w p hot equal to the velocity at maximum 
absorption in the Feii5169A line achieves high accuracy, 
as argued, e.g., by Schmutz et al. (1990). Unfortunately, 
this line is only present in synthetic spectra of models 
with relatively low T c g , corresponding to late times after 
explosion, when the use of either the SEAM or the EPM 
is not optimal fSoct. I2T3I4I . At early times, however, u p hoi 
can be inferred from hydrogen Balmer line measurements, 
but with a larger associated uncertainty that is difficult 
to estimate by eye-ball (Fig. 1141) . Thus, overall, the use of 
theoretical models like CMFGEN can contribute to reduc- 
ing the uncertainty on the inferred photospheric velocity 
by ca. 10-20% compared to methods based exclusively on 
direct measurement of observed line profiles, a significant 
gain given the desired accuracy for subsequent distance 
determinations. 

Despite the relatively good agreement obtained above 
(at the few tens of percent level), u a b s can, in some mod- 
els, be lower than w p hot- This is contrary to intuition. Since 
these lines have to be optically thicker than the continuum 
to leave an imprint on the emergent spectrum, one expects 
the line-flux minimum to occur exterior to, and thus at 
higher velocities than that at the photosphere. The pos- 
sible underestimate of i> p hot with the use of i> a bs can be 
inferred from the work of Branch (1977), and Jeffery & 
Branch (1990), but no explanation is given. 

To understand the behavior of the velocity minima, 
consider the example of the intermediate-stage model of 
Paper I (Sect. 3.2). In Fig. ED we see that the total op- 
tical depth of hydrogen Balmer lines (thick solid curves) 
is orders of magnitude larger than that in the continuum 
(solid or broken, thin curves). Yet, measurements on the 
synthetic spectrum (Fig. 3, Paper I) yield w a bs/w p hot of 
1.10, 0.99, and 0.93 for Ha, H/3, H7, confirming the typ- 
ical trend discussed in Fig. ^] We shall see below that 
the underestimate/overestimate of i> p hot is tied to the way 
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the emergent line intensity varies with depth and impact 
parameter p. 
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Fig. 16. Bottom: grayscale image of the quantity p ■ I(p) 
as a function of p and scaled wavelength x — (A/Ao — 
l)c/u p hot, where p is the impact parameter (in units of 
the hydrostatic radius R*) and I{p) the emergent specific 
intensity for a ray with p (at x). Here, Ao corresponds 
to the rest wavelength of Ha and c is the speed of light. 
The photospheric velocity of this model is 8750 km s -1 . 
The overplotted thick black curve gives the line-of-sight 
velocity location where the integrated continuum optical 
depth at 5500A, along z and at a given p, equals 2/3. 
Top: line profile flux, directly obtained by summing p-I{p) 
over the range of p. Overplotted are the curves p ■ I(p) for 
p = 0.05i?* and p = 1.8R* (thin lines; scaled arbitrarily 
and individually for clarity). 

Let us first focus on Ha. In the bottom panel of Fig. 1*11)1 
we show the gray scale image of the flux-like quantity 
p ■ I(p) as a function of scaled wavelength x = (A/Ao — 
l)c/w p hot and p. For clarity in this section, I(p) is com- 
puted by accounting only for the bound-bound transition 
in question. Overplotted (thick line) is a curve showing the 
location of the continuum photosphere, (£ P hot,pjP)- Along 
a given ray p, p ■ I{p) shows a minimum at x = x m i niP , 
with x m in !P < x p hot,p- This is to be expected because of 
the greater optical depth in the line than in the contin- 



uum, and because of the lower Ha source function (Paper 
I). However, following the curve (x p hot,p,p) as a guide, we 
find that the location x m i ntP varies greatly with p. This has 
important implications for the location of the absorption 
minimum in the observed profile. 

In the top panel of Fig. 1161 beside the total emergent 
line profile (solid line), we plot the quantity p- I(p) (scaled 
for clarity) for two selected impact parameters, p = 0.05i?* 
(impacting close to the SN center) and p = 1.8i?* (the lim- 
iting p where both line and continuum become optically- 
thin for all depths along the ray). For the inner ray, the 
minimum occurs at \x\ ~ 1.25, thus at a velocity well 
above w p hot; but for the outer ray, it occurs at \x\ ~ 0.8, 
thus for a velocity well below it. When building the total 
profile, one combines contributions from rays with impact 
parameters encompassing this range. The location of the 
minimum in the resultant profile will thus depend on the 
relative weighting of different impact parameters, which 
will be transition dependent. 

The cause of the variation in the location of the ab- 
sorption maximum was partially explained in Sect. 01 iso- 
velocity curves are at constant x (depth z in the (p, z) 
plane) but the density varies as l/r™. Thus, at fixed z, the 
density drops fast for increasing p, at the same time reduc- 
ing the probability of line scattering and/or absorption. In 
practice, the location of maximum absorption along a ray 
shifts to larger depths (z closer to zero) for increasing p, 
showing overall the same curvature as seen for the contin- 
uum photosphere (see also Fig. 1 1 0|> . 

In Figs. 1171181 we reproduce Fig. El f° r the cases of 
H/3 and Hy. The correspondence with Fig. El is striking, 
only mitigated by the difference in optical-depth (and line 
source functions) between these lines (Fig. IT5|) . Moving 
from Ha to H/3 and H7, the line optical depth decreases 
and the velocity at maximum absorption in the line profile 
shifts towards line center, i.e. from 1.1 to 1 and 0.95, in 
agreement with the values given above and measured on 
the synthetic spectrum (Fig. 3, Paper I), which includes 
all bound-bound transitions of all species. 

For completeness, using the late-stage model of Sect. 
3.3 (Paper I), we show in Fig. El the analogue of Fig. 1161 
for Fen 5169A (the oscillator strength of all bound-bound 
transitions of Fell located within 500A of the line cen- 
ter were set to zero). The maximum absorption occurs 
at \x\ ~ 1, and thus delivers exactly the photospheric 
velocity, in agreement with the velocity measurement of 
0.99w p hot obtained from the synthetic spectrum feature 
attributed to that line (Fig. 5, Paper I), which as above 
includes all bound-bound transitions of all species. 

Detailed model atmosphere calculations of type II 
photospheric-phase SN can thus deliver accurate and 
meaningful estimates of the photospheric velocity, by 
treating accurately the complex continuum and line op- 
tical depth effects controlling the line profile shape. The 
agreement between measurements done on line profiles as- 
suming either a single-line or all possible lines contribu- 
tions suggests that line-overlap does not affect noticeably 
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the inference of w p hot from the estimate of w a b s on the 
observed spectrum of photosphcric-phasc type II SN. 

A = 4340 A 




and p = 1.6.R*. Note that the velocity location of the flux 
minimum corresponds to a value smaller than w p hot- 

Fig. 17. Same as previous figure for H/3. Overplotted 
curves in the top panel correspond to rays with p = 0.05-R* 



and p = 1.8.R». Note that the velocity location of the flux umty . Compared to previous tabulations of correction fac- 
minimum coincides with u p hot- tors given by E96, our values are systematically higher, by 

ca. 0.1., although the reason is unclear. 

To gain further insight, we have extracted from these 

models the thcrmalization and photosphcric radii. The ra- 
6 Conclusion °^ these two radii is related to the dilution factor. We 

find that the ratio of the thcrmalization radius to the pho- 
In this paper, we have presented and discussed various as- tospheric radius is relatively uniform throughout the op- 
pects of the Expanding Photosphere Method, using knowl- tical and near-IR, and increases as the outflow cools down 
edge drawn from a large ensemble of CMFGEN models for and hydrogen recombines. For a given outflow ionization 
photospheric-phase type II SN. or effective temperature, the ratio is changed if the den- 

We have obtained such CMFGEN-based correction sity exponent or the scale of the supernovae is modified, 
factors from the minimisation technique at the heart of resulting from the modulation of the electron-scattering 
the EPM, used to compensate for the approximation of optical depth at the photosphere. We also find that if one 
the SN SED by a (unique-temperature) blackbody. The accounts for line opacity, the location of the photosphere 
correction factors do not just compensate for the effects radius can vary by up to a factor of 2-3, the effect be- 
of dilution due to electron-scattering - they also account ing maximum in the UV and more modest in the optical 
for the range of temperatures at the depth of thermal- and beyond. Nonetheless, this reveals that the uniqueness 
ization, and, most importantly, the presence of lines that of the photosphere radius for all wavelengths holds bet- 
alter considerably the SED at late times. This latter fact, ter at earlier times, when type II SN spectra are more 
usually ignored by other groups, is what causes correction featureless, closer to the equivalent continuum energy dis- 
factors to take values of up to 1.5 in our grid of models, tribution. This has some bearing on the adequacy of using 
while the upper limit for dilution factors is, by definition, the Baade method for distance determinations. 
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X = 5169 A 




(A/A -1)c/v 

* ' ' ' phot 

Fig. 19. Same as previous figure for Feil5169A, in the cool 
model presented in Sect. 3.3 of Paper I, corresponding to 
a late-stage model in which Fe II lines are present through- 
out the V band. For this plot, we reset the oscillator 
strength of all other Fell contributors to this wavelength 
region to essentially zero, in order to highlight the forma- 
tion region of the single transition at 5169A. Overplotted 
curves in the top panel correspond to rays with p = 0.04i?„ 
and p = 1.5.R*. Note that the velocity location of the flux 
minimum corresponds to a value consistent with u p hot, just 
like for H/3 in Fig. El f° r an intermediate-stage (hotter) 
model. 



In addition to the determination of correction factors, 
the accuracy of the measurement of the photospheric ve- 
locity is equally important for the reliability of the EPM. 
We have shown that, contrary to expectations, optically 
thick lines do not necessarily show a minimum P-Cygni 
line profile flux at a Doppler-shifted wavelength that cor- 
responds to the photospheric velocity. Instead, depending 
on the outflow properties, such a measurement can deliver 
an overestimate or an underestimate of the photospheric 
velocity. This is particularly problematic for earlier mod- 
els which show broad P-Cygni line profile troughs, mostly 
for hydrogen Balmer lines. Unfortunately we have also 
demonstrated that, due to the more well defined photo- 
spheric radius, the lack of contaminating lines and a SED 



closer to that of a blackbody, it is at these earlier times 
that the EPM is best used. 

This paper has established some foundations for a care- 
ful use of the EPM, which we are currently applying to 
SN1999em. 
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